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ABSTRACT 

We propose that luminous transients, including novae and supernovae, can be used to detect the faintest 
galaxies in the universe. Beyond a few Mpc, dwarf galaxies with stellar masses < 10® Mq will likely be 
too faint and/or too low in surface brightness to be directly detected in upcoming large area ground-based 
photometric surveys. However, single epoch LSST photometry will be able to detect novae to distances of 
~ 30 Mpc and SNe to Gpc-scale distances. Depending on the form of the stellar mass-halo mass relation and 
the underlying star formation histories of low mass dwarfs, the expected nova rates will be a few to ~ 100 yr“' 
and the expected SN rates (including both type la and core-collapse) will be ~ 10^- 10"* within the observable 
(Ttt sr) volume. The transient rate associated with intrahalo stars will be comparably large, but these transients 
will be located close to bright galaxies, in contrast to the dwarfs, which should trace the underlying large 
scale structure of the cosmic web. Aggressive follow-up of hostless transients has the potential to uncover the 
predicted enormous population of low mass field dwarf galaxies. 

Subject headings: galaxies; stellar content — galaxies; dwarf — supernovae; general — novae, cataclysmic 
variables 


1. INTRODUCTION 

Dwarf galaxies are believed to play several important roles 
in modern theories of galaxy formation. The ionizing radi¬ 
ation they emit at high redshift is likely es s ential for reion¬ 
izing the universe (e.g., iMadau et al.lll999l; iRobertson et al.l 
120131) . At later times they contribute to the buildup of 
streams and stellar hal os around galaxies (e.g.. lJohnstonll99^ 
iHelmi & WhitellT999h . Their luminosity functions and kine¬ 
matic properties place very demanding constraints on mod¬ 
els of star formation and stellar feedback and even perhaps 
the nature of dark matter (e.g.. iBovlan-Kolchin et al.l 1201 It 

iBrooks et al.ll2013HGovernato et al.ll2015l). 

The relation between stellar mass and dark matter halo 
mass is a simple but powerful tool to help understand the 
distribution of galaxies in the cosmic web, the luminos¬ 
ity function of dwarf galaxies, the integrated efficiency of 
star formation, and the hierarchical growth of galaxies over 
cosrn i c time (e.g.. iPurcell et al.l l2007t IConrov & Wechslej 
l2009t iBehroozi et akl ^13t iGarrison-Kimmel et al.l l2014l) . 
Empirically-constrained stellar mass-halo mass relations im¬ 
ply that there exist enormous numbers of dwarf galaxies 
with stellar masses < 10®in the field (> 1 Mp c~^; e.g., 
IBehroozi et al.ll2013l; iGarrison-Kimmel et al.l 120141) . Only a 
handful of these v ery low mass field dwarf s have been de¬ 
tected to date (e.g.. lKarachentsev et al.lUOl^ . 

Their low luminosities and surface brightnesses makes 
dwarf galaxies difficult to detect. Dwarf galaxies within a few 
Mpc can be resolved into individual stars. This has been the 
class ic detection method for nearby dvyarfs in the modern era 

(e.g JWi llman et al.l2005HMartin et al.l2006l;lBelokuroy et al.l 

l2007t Ichiboucas et al.lboool) . The faintest dwarfs haye sur- 

face brightness within a half light radius of > 27 mag arcsec”^ 
(lMcConnachiell20T^ . making it yery difficult to detect them 
at distances where they cannot be resolyed into stars. More- 
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oyer, there are theoretical reasons to suspect that a large popu¬ 
lation of ultra-diffuse galaxies may exi st beyond the detection 
limits of current res olyed-star suryeys (iBoyill & Ricottil[2009l 
iBullock et al.ll2010i Wheeler et ak, in preparation). Special¬ 
ized imagers can reach these faint surface brightness l imits 
(e.g..lMih ^os et al.l l2005l;lAbraham & yan Dokkumll2014h . and 

iMerritt et al.l (120141) recently reported the discoyery of faint 

dwarf galaxies surrounding MlOl detected in this way. How- 
eyer, detecting faint dwarfs in the field requires imaging large 
areas of the sky, and no such suryeys targeting these extremely 
low surface brightness limits are currently planned, though 
clearly they would be of high yalue. At still greater distances 
dwarfs will appear as faint point sources, and so either photo¬ 
metric redshifts or follow-up spectroscopy is required to con¬ 
firm their dwarf status. 

Ongoing and planned high cadence photometr ic suryeys 
such as the Palomar Transient Eac tory (PTE; Rau et akl 
[200l;lLawet'ani200^ . Pan-ST ARRS (iKaiser et akll2002^ the 
Zwicky Transient Eactory (ZTE iBellrn 201A ), and the Large 
Synoptic Suryey Telescope (LSST; lAbell et alll2009h are or 
will be obserying large areas of the sky with the depth neces¬ 
sary to detect a wide class of transient objects in single epoch 
imaging and indiyidual bright stars in nearby dwarf galaxies 
in coadded images. In this Letter we propose that the tran¬ 
sients that will be discoyered in wide field photometric sur¬ 
yeys can be used to detect faint dwarf galaxies that would 
otherwise be too distant to detect with standard techniques. 
This proposal is similar in spirit to the idea of using transients 
to detect the diffuse population of stars in the halos of galax- 
ies, groups, and clusters, often known as intrahalo st ar s (e.g., 
iGal-Yam et akl 120031 [Sh^lT?)?)^ Is^et akl 1201 2h . [Sh^ 
(I2006h for example predicted that LSST will obserye hun¬ 
dreds of noyae associated with intrahalo stars. 


2. DWARF GALAXIES IN THE LOCAL GROUP 

iMcConnachiel (1201 2h proyided a compilation of key photo¬ 
metric and structural paramaters of all known dwarf galaxies 
(as of 2012) within 3 Mpc of the Sun. Erom the reported ef- 
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Figure 1. Properties o f know n dwarfs in the Local Group, from data com¬ 
piled bv lMcConnachid 42012f) . The dotted line in the upper panel marks the 
rough expected surface brightness limit of LSST. Stellar masses have been 
estimated assuming M/Ly = 2.0. 


fective radii, magnitudes, and ellipticities we have computed 
the mean V-band surface brightness within the effective ra¬ 
dius, /Teff v- The distribution of known Local Group dwarf 
galaxies in size-luminosity-surface brightness space is shown 
in Figure[T] We also translate the sizes into distances at which 
the size would subtend 1". A threshold surface brightness 
of f^eff.v = 27 mag arcsec“^ is indica ted as a rough est imate 
of the LSST surface brightness limit (lAbell et alJl2009l) . No¬ 
tice that the faintest known galaxies all lie at ^eff.v > 27 mag 
arcsec“^, suggesting that it will be difficult for LSST to de¬ 
tect such galaxies except where the individual stars from such 
galaxies can be easily detected (e.g., within a few Mpc of the 
Sun). In the figure we also translate the x-axis into an ap¬ 
proximate stellar mass by assuming a V-band mass-to-light 
ratio of M/Ly = 2.0, appropriate for old (10 Gyr) metal-poor 
([Z/H] < -1) populations (IConrov et al.l2009Pl Of course, for 
the star-forming galaxies this is a poor assumption and will re¬ 
sult in an overestimation of the mass. These masses are only 
computed to guide the eye and are not used in any of the cal¬ 
culations below. 

^ A IKroup^ 120011) stellar initial mass function has been assumed through¬ 
out this work when quoting stellar masses and rates per unit stellar mass. 


log n(>M|,) Mpc ^ 

3.41 1.83 0.00 -1.5 -2.6 



Figure 2. Stellar mass-halo mass relations. The upper x-axis has been con¬ 
verted to cumulative number density based on our fiducial cosmology. The 
B13 relation is a significant extrapolation beyond the mass range over which 
the original relation was constrained. The shaded region marks the approx¬ 
imate region where reionization is expected to quench star formation in low 
mass halos. 


3. MODEL INGREDIENTS 

In this section we present a model for the expected SN and 
nova rates from dwarf galaxies and intrahalo stars. 

We start with the z = 0 dark matter halo mass function from 
[Tinker et alJ (l200 8h. The halo virial mass is defined accord- 
ing to lBrvan & NormanI (Il998h . We adopt a Hubble constant 
of Hq = 70 km s“' Mpc“\ U„, = 0.27 = 1 - Ha, and q-c = 0.8 , 
consistent with the 7-yr WMAP results (iKomatsu et alj20Tl]l . 
These parameters are somewhat different from the r ecent 
Planck measurements dPlanck Collaboration et al.ll2015l) . but 
the differences have no material impact on our results. 

Galaxies are placed within dark matter halos according to 
a stellar mass-halo m ass (M^-Mi,) relation. Our fid ucial re¬ 
lation is adopted from iGarrison-Kimmel et aP (1201 4 GK14), 
which reproduces both the empirical galaxy stellar mass func¬ 
tion at z = 0 and the Local Group dwarf luminosity function. 
The form of the M* - M/, relation at low masses is uncer¬ 
tain and has a large effect on our model predictions. In or¬ 
der to expl ore this sensitiv i ty we will also consider the re¬ 
lation from iBehroozi et ^ (120131 B13), which we extrap¬ 
olate well beyond the mass range probed by B13 (roughly 
M* = 10^- 10 ^'Mq). The two relations are compared in Fig¬ 
ure 121 We also indicate with a shaded band the approximate 
region where reionization i s expected to quenc h star forma¬ 
tion in low mass halos ('e.g.. lBullock et al1l2000h . We adopt a 
relation with no scatter whereas i n reality there may in fact be 
significant scatter at low masses (I Wang et alj|2015l Garrison- 
Kimmel et al., in preparation). As long as the scatter is sym¬ 
metric about the mean relation this will not affect our pre¬ 
dicted rates. 

The next step is to specify star formation histories (SFHs) 
for the dwarfs, as the SFHs will influence the SN rates. 
We adopt the following simple model: dwarf galaxies with 
M* > 10^ Mq are assumed to have a constant SFR through¬ 
out their lifetimes, while less massive dwarfs are assumed to 
form all of their stars instantly at high redshift (z > 2). We 
consider in this work only dwarfs in the field, i.e., dwarfs that 
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Figure 3. Predicted SN and nova rates as a function of distance adopting the GK14 (left) and B13 (right) stellar mass-halo mass relations. Lines show rates for 
dwarf galaxies with stellar masses less than IO^Mq, IO^Mq, and 10^ Mq. The shaded regions indicate where SNe and novae are detectable with a single LSST 
visit (r < 24.3) and where the rates are > 1 yr“^. Each line is broken into three regions: where a typical galaxy would be spatially resolved (dotted), unresolved 
and brighter than r = 27 (solid) and both unresolved and fainter than r = 27 (dashed). This depth was chosen to be comparable to the final 10 yr LSST coadded 
depth of 27.5. Note that the SN/nova ratio is different for the lowest mass bin compared to the higher mass bins owing to the different assumed SFHs. The 
quoted nova rates on the y-axes ai‘e therefore only appropriate for the two most massive bins. 


are not within the virial radius of a larger halo. This model 
is motivated by SFHs derived from res olved color magni tude 
diagrams of Local Group dwarfs (e.g., IWeisz et^l2014ll and 
also consistent with expectations for dwarfs above and below 
this stellar mass thr eshold in cosmological simulations (e.g., 
lOnorbe et al.l 1201 M Wheeler et al. in preparation). Later, 
when converting masses into r-band magnitudes, we adopt 
M/Lr = 1.8 for the low mass galaxies with ancient metal- 
poor stellar populations and M/Lr = 0.73 for th e higher mass 
dwarfs with constant SFHs (IConrov et al.ll2009h . 

We will genetically refer to SNe, noting that this includes 
both the core-collapse (CC) and type la varieties. For the CC 
SN rate we assume 1 per IOOMq of newly formed stars, i.e., 
a galaxy with a star formation rate of 1 Mq yr“^ will have CC 
SN rate of 10“^ yr“'. It is well-known that the rate of la SN 
depends on the age of the progenitor stellar population, with 
a delay time distribution (DTD) that falls approximately as 
f~*. W e adopt a DTD based on data presented in iMaoz et alJ 
([20T2). with a rate at 10 Gyr of 3 x 10“^"^ yr * Mq and a time 
dependence of . With this DTD, a galaxy with a constant 
SFH will at the present epoch have a la rate of « 10“'^ yr“' 
Mq (assuming the DTD is zero at f < 10* yr). The above 
assumptions imply that the CC rate is lOx higher than the la 
rate for a galaxy with a constant SFH. 

Distances to the SNe are required to use them to detect 
low redshift dwarf galaxies. For type la SNe, deriving dis¬ 
tances from the light curve is straightforward. There has been 
some effort to use a su bset of CC SNe (the “plateau” type) as 
stand a rd candles (e.g., ISchmidt et all ll994tlHamuv & Pinto! 
120021) . iHamuv & Pintol (120021) showa reasonably strong cor¬ 
relation between type IIP SNe apparent magnitude and red- 
shift, which should be sufficient to at least discriminate be¬ 
tween high and low redshift SNe. Spectroscopic follow-up 
of transients, which will no doubt occur for a subset of the 
events, will also be valuable for confirming the redshift of 
these apparently hostless SNe and therefore also the redshift 
of the underlying faint dwarf galaxy. 

Nova rates have historically been difficult to empirically 
estimate owing to a variety of selection and incomplete¬ 


ness effects. Nonetheless, observations of galaxies in the 
nearby universe are broadly consistent with a universal nova 
rate of 2 x 10~*° yr~^ML*, independent of Hubble type , 
galaxy color, etc. (ICiardullo et alJl990HFerrarese et all2003l : 

IWilliams & Shafteill2004l) . There is however some evidence 

that the nova rate m ay increase toward lower mass galaxies 
(iNeill & Sharall2005h . which, if confirmed, would haye im¬ 
portant implications for the present study. Noyae peak abso¬ 
lute magnitudes fall in the range from -7 < My < -9; in this 
work we adopt a typical peak magnitude of Mr = -8. The 
decay time of the noya light curye is inyersely related to its 
peak brightness and the relation is sufficiently tight, and the¬ 
oretically well-understood, so that one may be able to infer 
distances to noyae from th eir light curyes (e.g.. lSharalll981l : 
iDownes & Duerbe^l2000h . 

iPurcell et alJ (l2007b~ presented a simple model for the 
buildup of intrahalo light (IHL), where intrahalo stars are de¬ 
fined as residing within a parent dark matter halo but not 
bound to any particular galaxy. The model assigns stars to 
dark matter halos in a cosmological N-body simulation yia 
an empirically-constrained M* -M/, relation, and when a dark 
matter subhalo is disrupted in the simulation the stars from 
that subhalo are assigned to the IHL. At the cluster scale, 
this model produces IHL fractions of order unity while at the 
Milky Way halo scale the IHL contributions are at the percent 
leyel; both of these predictions are in agreement wi th obserya- 
tions re.g.. lGonzalez et ^l2005Hlrwin et alJl2005l) . We haye 
fit their results with a polynomial function: 

M ^ 

/,HL(M„) = -;7^ = ^fl/(logM,)', (1) 

with a, = (9.5936,-5.3943,0.5965,-0.01869), where M*,ihl 
and M* gal are the stellar masses in the IHL and the central 
galaxy, respectiyely. The fit reproduces the model results to 
5- 10%. We assume that the intrahalo stars are uniformly 
old when computing noya and SN rates from this component. 
Integrating oyer the halo mass function, the aboye equation 
implies that the total stellar mass in the IHL is Ri 10% of the 
















































4 


CONROY & BULLOCK 


stellar mass in galaxies. 


4. RESULTS 

The resulting model SN and nova rates for dwarf galaxies 
are shown in Figure[3as a function of distance. We show re¬ 
sults for two stellar mass-halo mass relations in order to high¬ 
light the sensitivity to the adopted relation. Clearly, observa¬ 
tions of the SN and nova rate could place powerful constraints 
on the low mass end of the M* -M/, relation. Rates are shown 
for dwarfs with stellar masses less than IQ^Mq, 10® Mq, and 
10^Mq. Along each line we mark where the most mas¬ 
sive dwarfs would be resolvecQ, unresolved and brighter than 
r = 27, and unresolved and fainter than r = 27. This depth was 
chosen to be comp arable to the fina l, 10 yr coadded LSST 
depth of r 27.5 (lAbell et alJl2009t) . In order to compute 
when a dwarf would be resolved we have assumed effective 
radii of 0.8, 0.3, and 0.2 kpc for dwarfs with stellar masses 
of 10*^ M 0 , 10® Mq, and 10® M 0 , respectively (see Fig [T]i. 
Shaded regions denote “sweet spots” where the rates are > 1 
yr“^ and the transients are observable. In the case of novae 
the maximum distance is set by being able to detect a = -8 
nova in single epoch LSST data (though of course multiple 
epochs will be required to confirm tha t it is a nova), wh ich 
will have a minimum depth of r Ri 24.3 (lAbell et alJl2009h . In 
the case of SNe the maximum depth is set by being able to 
detect a Mr = -18 SN in single epoch LSST data. This is a 
conservative limit as both type la and CC SNe can have peak 
magnitudes brighter than this. Finally, note that the quoted 
rates are for the entire (Ttt sr) observable volume, whereas in 
practice any survey will cover only a fraction of the sky at the 
necessary cadence. 

The principle result of this Letter is that SNe and novae 
residing in low mass dwarf galaxies should be detectable in 
large numbers with upcoming time domain surveys such as 
LSST. Dwarfs in the mass range of M* ~ 10® -10®M 0 should 
be spatially resolvable to 10-100 Mpc distances, but their 
low surface brightness, fi > 26 mag arcsec“^, will make it 
challenging to directly detect them. However, novae within 
these low mass galaxies will be detected at rates of 1-10^ 
yr“* within ^ 30 Mpc, and SNe will be detected out to Gpc 
scales at rates of ^ 10^-10^ yr“* (note that a comoving radial 
distance of 1 Gpc corresponds to z = 0.25 in our adopted cos¬ 
mology). Even very low mass dwarfs with M* < 10® M 0 will 
produce type la SNe at a rate of ~ 10^ yr“*. Relatively mas¬ 
sive dwarfs (M* ^ 10® M 0 ), should be easily detectable as 
spatially resolved, relatively high surface brightness objects 
out to several hundred Mpc distances. For these objects, de¬ 
tection via transients will likely be less interesting. 

Figure |4] shows the expected SN and nova rates from in¬ 
trahalo stars. Results are shown for stars within dark matter 
halos less massive than 1O'®M0 and 1O^^M0. The implied 
nova rates within 30 Mpc are fairly high, reaching nearly 
« 500 yr“* if one considers all IHL stars. A key difference 
between the transients associated with IHL stars and dwarfs 
is their spatial distribution - IHL stars are confined to within 
0.1-1 Mpc (depending on the virial radius of the host halo), 
while the dwarfs are presumably distributed throughout space, 
tracing the underlying dark matter distribution with a con- 

Here we use “resolved” and “unresolved” to refer to spatially resolv- 
ing the galaxy, as opposed to resolving the galaxy into individual stars. We 
consider a galaxy “resolved” if the effective radius is > 0''.5. 



10 ** 



10 - 


10 " 


Figure 4. Same as Figure[3] except now for intrahalo stars residing in halos 
less massive than 10*^ Mq and 10*^ Mq. The intrahalo stars are assumed to 
be uniformly ancient and so the SN/nova ratio is different here than in the 
y—axes in Figure[^ 

stant anti-bias. The implied IHL rates are broadly consis- 
tent w i th, alt hough somewhat lower than the predictions from 
ISharal (l2006h . Shara assumed a global IHL fraction, /ihl, of 
0.1 -0.4, whereas our model for the IHL, constrained by re¬ 
cent data, places this fraction at closer to 0.1. 

5. DISCUSSION 

Dwarfs at the mass scale M* ^ 10® - 10®M 0 reside within 
some of the smallest dark matter halos that are expected to 
form stars - scales where the physics of reionization and dark 
matter may be probed with potentially powerful effect. Un¬ 
fortunately, their use as cosmic laboratories to this point has 
been fundamentally limited by the fact that we detect them 
with reasonable completeness only within ~Mpc volumes. 
As Figure[3]makes clear, the detection of (apparently) hostless 
SNe at 0.1 - 2 Gpc distances should provide an important av¬ 
enue for detecting these dwarfs and constraining the M* -M/, 
relation at low masses within volumes that are virtually im¬ 
mune to sample variance effects. In addition to counts, the 
ratio of type la to CC SNe will be a sensitive probe of the av¬ 
erage SFH of the unresolved dwarfs. The correlation function 
of hostless SNe within a Gpc will provide constraints on the 
halo masses of the dwarfs and will allow for a statistical sep¬ 
aration of transients arising from dwarfs and from intrahalo 
stars. One could also imagine analyzing the rates as a func¬ 
tion of large scale environment to probe the low mass end of 
the galaxy population in e.g., void regions. 

The majority of the dwarf galaxies detected via transients 
will be too distant for detailed study. We envision this detec¬ 
tion technique as primarily a statistical probe of the low mass 
dwarf population. A dwarf with M* = 10® M 0 and a constant 
SFH will have a CC SN rate of 10“® yr“', so one must observe 
a volume containing 10® such dwarfs to detect one SN per 
year. This technique is only feasible because of the enormous 
numbers of predicted dwarf galaxies in the field. 

A number of assumptions were made in computing the rate 
of luminous transients, the most important of which is the 
adopted relation. Other important assumptions in¬ 

clude the SFH of the dwarfs and the assumed constant nova 
rate per unit stellar mass. We have assumed that galaxies 
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with M* > 10^ Mq have constant SFHs in the field. If in¬ 
stead some fraction are quenched or the star formation rates 
were higher in the past then the predicted SN rates will be 
lower. On the other hand, if the nova rate is in fact h igher 
in dwarf galaxies, as suggested bv iNeill & Sharal (120051) . then 
the expected nova rate from dwarf galaxies could be much 
higher than assumed here. LSST and other transient surveys 
will directly constrain the nova rate in nearby dwarf galax¬ 
ies with well-studied stellar populations, which will effec¬ 
tively eliminate the nova rate as a source of uncertainty. We 
have also assumed that the CC SN rate per unit star forma¬ 
tion rate is constant and does not change with metallicity nor 
with the level of star formation rate. The latter in particular 
might affect the SN rate at levels where stochasticity and/or 
the effects of small cloud masses become im portant (e.g., 
IWeidner & Kroupal 120041 : [da Silva et al.ll201^ . Here again 
LSST and other transient surveys will provide rates of CC 
SNe within known galaxies to inform a better model for the 
expected CC SN rate in low mass dwarfs. Finally, we have 
included CC SNe in the transient budget, but it is not clear 
if these will be as useful as novae and type la SNe due to 
the difficulty in assigning distances to CC SNe based only 
on the light curve. However, even in the most pessimistic 
case in which one discards the CC SN events, the predicted 
la rates forM* > 10^ Mq are still substantial (approximately 
10 X lower than shown in Figure |2l. Moreover, recall that our 
model for the lowest mass galaxies (M* < IQ^Mq) assumes 
an ancient stellar population and hence the SNe associated 
with those galaxies are solely of the la variety. 

We thank Jeremy Tinker for providing code to compute halo 
mass functions and Mike Boylan-Kolchin, Mansi Kasliwal, 
Pieter van Dokkum, and Dan Weisz for comments on an ear¬ 
lier draft. CC is supported by a Packard Foundation Fellow¬ 
ship. This research was supported in part by the National Sci¬ 
ence Foundation under Grant No. NSF PHYl 1-25915. The 
authors acknowledge the KITP workshop "Galactic Archae¬ 
ology and Precision Stellar Astrophysics" where this work 
was conceive d. Ned Wright’s online cosmology calculator 
(IWrigh)i2006l) was used in the preparation of this manuscript. 
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